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Abstract 

Resummation of large QCD radiative corrections, including leading and 
next-to-leading logarithms, in pion electromagnetic form factor is reviewed. 
Similar formalism is applied to exclusive processes involving heavy mesons, 
and leads to Sudakov suppression for the semi-leptonic decay B -^ nlu. It is 
found that, with the inclusion of Sudakov effects, perturbative QCD analysis 
of this decay is possible for the energy fraction of the pion above 0.3. By 
combining predictions from the soft pion theorems, we estimate that the 
upper limit of the KM matrix element \Vub\ is roughly 0.003. 



1. Introduction 

It has been shown that perturbative QCD (PQCD) is apphcable to ex- 
clusive processes such as elastic hadron form factors for energy scale higher 
than few GeV [|l|. The enlargement of the range of the applicability from 
much higher energies [@] down to this low scale is due to the inclusion of 
transverse momentum dependence into factorization theorems. This depen- 
dence appears in an exponential factor, which arises from the resummation 
of large radiative corrections, and leads to Sudakov suppression for elastic 
scattering of isolated colored quarks. Detailed derivation of Sudakov fac- 
tors for hadron-hadron Landshoff scattering refers to [Q]. Similar expressions 
have been obtained and employed in the PQCD analysis of pion and proton 
form factors |l|, H, |^, H , pion Compton scattering [0 and other exclusive pro- 
cesses [^. Resulting predictions from this modified version of factorization 
theorems have been examined and found to be dominated by perturbative 
contributions. 

All the above analyses of large corrections involves only light hadrons. In 
this paper we shall extend the formalism to exclusive processes containing 
both light and heavy mesons, such as B meson decays, and organize their 
Sudakov corrections to all orders. An important work in the study of the 
standard model is to determine the mixing angles in the Cabibbo-Kobayashi- 
Maskawa (KM) matrix. The decay K — > ttIu contains the information of the 
matrix element \Vus\, and chiral symmetry provides a precise method to study 
this process ||^. Similarly, \Vcb\ is determined by exploring the B -^ Dlu 
decay, for which heavy quark symmetry is an appropriate tool |Ty]. As to 
\Vub\i which can be measured reliably from the B -^ ttIu decays ||lT|, neither 
of the above theories is proper. 

Recently, an analysis of the semi-leptonic decay B —>■ nlu based on the 
heavy quark effective theory (HQET) has been performed by Burdman et al. 
[O . They have determined the normalization of form factors involved in the 



process in terms of the soft pion relations, and given the ratio of these form 
factors to the corresponding ones in the D -^ ttIu decay. However, exphcit 
evaluation of this process is not yet successful. In this paper we shall show 
that, by incorporating Sudakov effects, PQCD is indeed applicable to the 
semi-leptonic B meson decays as the pion carries away more than a quarter 
of its maximum energy. The differential decay rate is then obtained. By 
combining our predictions with those from the soft pion theorems, which 
have been derived in the framework of HQET [|12|, the total decay rate and 
the branching ratio of this process are estimated. Comparing this estimation 
with experimental data, we find that the upper limit of the KM matrix 
element \Vub\ is roughly 0.003. PQCD is then able to complement HQET in 
the analysis of heavy meson decays. Our results can be easily generalized to 
perturbative analysis of other heavy-to-light transitions. 

We find that factorization theorems are successful for this semi-leptonic 
decay especially when the pion is energetic. In this case non-perturbative 
approaches such as the soft pion theorems and QCD sum rules are not appli- 
cable. The first attempt to apply the modified perturbative formalism to the 
analysis of the B -^ ttIu decay has been made by Akhoury et al. [|13|. How- 
ever, they did not consider the transverse momentum dependence, and have 
obtained predictions which are drastically different from ours. As claimed in 
[ |l^ , the hard gluon involved in the decay process, with Sudakov effects taken 
into account, is off-shell at most by an amount of 1.4AQCD'^fe; ""^fe being the h 
quark mass. In our approach the hard gluon is off-shell roughly by SAqcD^^fe, 
and therefore, the perturbative analysis is more reliable. 

The resummation of Sudakov logarithms for the pion form factor is re- 
viewed in section 2. The same formalism is then applied to semi-leptonic B 
meson decays in section 3, where the full expression for the Sudakov factor 
including leading and next-to-leading logarithms is given. Section 4 contains 
numerical analysis of the modified factorization formulas for the relevant 
form factors and the differential decay rate oi B ^ -kIv. Section 5 is our 



conclusion. 

2. The Pion Form Factor 

First, we review the modified factorization formula for a simple light-to- 
liglit process, the pion electromagnetic form factor, to lowest order of cou- 
pling constant a^, which is expressed as the convolution of wave functions 
and a hard scattering amplitude. We then investigate radiative corrections 
to this formula, and explain how they are absorbed into the above convolu- 
tion factors. The first step is to find out the leading momentum regions of 
radiative corrections, from which important contributions to loop integrals 
arise. There are two types of important contributions: coUinear, when the 
loop momentum is parallel to the incoming or outgoing pion momentum, 
and soft, when the loop momentum is much smaller than the momentum 
transfer Q^ of the process. Here, Q^ is assumed to be large and serves as an 
ultraviolate cutoff of loop integrals. We associate small amount of transverse 
momenta k^ with the partons in the above factorization picture, which is 
taken as an infrared cutoff. 

Each type of these important contributions gives a large single logarithm. 
They may overlap to give a double (leading) logarithm in some cases. These 
large logarithms, appearing in a product with a^, must be organized in order 
not to spoil the perturbative expansion. It is known that single logarithms 
can be summed to all orders using renormalization group methods, and dou- 
ble logarithms must be organized by the technique developed in ||I4[. We 
shall work in axial gauge n ■ A = 0, in which the resummation technique is 
developed easily, n being the gauge vector and A the gauge field. 

The diagrams shown in fig. 1 represent the 0[a^ radiative corrections to 
the basic factorization of the pion form factor, which contain large logarithms 
mentioned above. In axial gauge the two-particle reducible diagrams, like 



figs, la and lb, have double logarithms from the overlap of collinear and soft 
enhancements, while the two-particle irreducible corrections, like figs. Ic and 
Id, contain only single soft logarithms. This distinction is consistent with the 
physical picture: two partons moving in the same direction can interact with 
each other through collinear or soft gluons, while those moving apart from 
each other can interact only through soft gluons. Below we shall concentrate 
on reducible corrections, and demonstrate how they are summed into the 
Sudakov factor. 

A careful analysis shows that soft enhancements cancel between figs, la 
and lb, as well as between Ic and Id, in the asymptotic region with b ^ 0, 
b being the conjugate variable to ky. Therefore, reducible corrections are 
dominated by collinear enhancements, and can be absorbed into the pion 
wave funtion, which involves similar dynamics. Irreducible corrections, due 
to the cancellation of their soft divergences, are then absorbed into the hard 
scattering amplitude. Hence, the factorization picture still holds at least 
asymptotically after radiative corrections are included. The above cancella- 
tion of soft divergences is closely related to the universality of wave functions. 
For large b, double logarithms are present and the resummation technique 
must be implemented. 

Based on the above reasoning, the factorization formula for the pion form 
factor in the b space is written as |jl|, 

/■I r d b 

KiQ'^) = J dxidx2 I -——^V{x2,b,P2,fi) 

X H{xi,X2,b,Q,fi)V{xi,b,Pi,fi) , (1) 

where fi is the factorization and renormalization scale, and b is the separation 
between the two valence quarks. Pi and P2 are momenta of the incoming 
and outgoing pions, respectively. Here we choose the Breit frame, in which 
P^ = P2 = Q/V^ and all other components of P's vanish, Q^ being the 
momentum transfer, Q^ = —(Pi — P2)^. Note that eq. (^ depends only on a 



single b, because the virtual quark line involved in the hard scattering H is 
thought of as far from mass shell, and is shrunk to a point |1| . The pion wave 
function V includes all leading logarithmic enhancements at large b. The 
basic idea of the resummation technique is as follows. If the double logarithms 
appear in an exponential form V ~ exp[— const, x In Q ln(ln Q/ In b)], the task 
will be simplified by studying the derivative of V, dP/d InQ = CV. It is 
obvious that the coefficient C contains only large single logarithms, and can 
be treated by renormalization group methods. Therefore, working with C 
one reduces the double-logarithm problem to a single-logarithm problem. 

The two invariants appearing in V are P ■ n and n^, and by the scale 
invariance of n in the gluon propagator, 

M^^(,^ = =lL.^-"^£±^^n'J^], (2) 

q \ n ■ q [n ■ q) J 

V can only depend on a single large scale z/^ = (P ■ nf' jr? . It is then easy 
to show that the differential operator d/d In Q can be replaced by d/dn: 

^ -P = -^P°Ap. (3) 



lO 



dlnQ P-n dn° 

The motivation for this replacement is that the momentum P flows through 
both quark and gluon lines, but n appears only on gluon lines. The analysis 
then becomes easier by studying the n, instead of P, dependence. 
Applying d/drio, to the gluon propagator, we get 

^ -AT/^^ = L(Ar/^«g'^ + N'^'^q^) . (4) 



dn^ q ■ n 

The momentum q that appears at both ends of a gluon line is contracted with 
the vertex, where the gluon attaches. After adding together all diagrams with 
different differentiated gluon lines and using the Ward identity, we arrive at 
the differential equation of V as shown in fig. 2a, in which the square vertex 
represents 

P ■ nq ■ n 
6 



T" being the color matrix. An important feature of the square vertex is that 
the gluon momentum q does not lead to coUinear enhancements because of 
the nonvanishing n^. The leading regions of q are then soft and ultraviolet, 
in which fig. 2a can be factorized according to fig. 2b at lowest order. The 
part on the left-hand side of the dashed line is exactly V, and that on the 
right-hand side is assigned to the coefficient C introduced before. 

Therefore, we need a function /C to organize the soft enhancements from 
the first two diagrams in fig. 2b, and Q for the ultraviolet divergences from the 
other two diagrams. The soft substraction employed in Q is to avoid double 
counting. Generalizing the above two functions to all orders, we derive the 
differential equation of P, 



P(x,6,P,/i) 



dlnQ 



2/C(6/i) + ^ g{xvlii) + 1 ^((1 - x)u/ii) 
xV{x,b,P,fi). (5) 



The functions /C and Q have been calculated to one loop, and the single log- 
arithms have been organized to give their evolutions in b and Q, respectively 
p. They possess individual ultraviolet poles, but their sum /C + Q/2 is finite 
such that Sudakov logarithms are renormalization-group invariant. 

Substituting the expressions for /C and Q into eq. (j^), we obtain the 
solution 



V{x, b, P, /u) = exp 



P(a;,6,/i). (6) 



^=x, 1—x 

The exponent s, grouping the double logarithms in V, is expressed in terms 
of the variables 



asO 



q ^ ln[eQ/(v^A)] 

b = ln(l/6A) (7) 

A«.. fq\ A(2) fq ^ A« 



.«,,g) . _„„|U^ 1-1 -_(,-. 



A(i)/?2 . rin(26) + 1 

Tepf'^i b 

4/3? 4A ""l 2 



ln(2g) + 1' 


\ ■ 


V6/ 



32/5? 



In2(2g) - ln2(2S) 



(8) 



with A = Aqcd- The coefficients A^^^ and /3j are 

33 - 2n; 153 -19n; 

Pi — — tt; — 5 P2 — - 
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^(2) 



24 
67 TT^ 10 



-n 
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-A In , 
3 V2 



(9) 
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where tt,/ = 3 in this case is the number of quark flavors and 7 is the Euler 
constant. To derive eq. (H), the choice of the gauge vector n'^ oc (Pi + P2)'^ 
has been made. 

The function V{x, b, fi) and H still contain single logarithms from ultravi- 
olet divergences, which need to be summed using their renormalization group 
equations B: 



VVix,b,fx) = -2-fgV{x,b,fx) 

V H{xi, b, Q,fi) = 4 7g H{xi, b, Q, fi) 



(10) 



with 






(11) 



7g = —as/71 is the quark anomalous dimension in axial gauge. Solving 
eq. (PUJ), the large-6 asymptotic behavior of V can be summarized as 



V{x, 6, P, fj) = exp 



xV{x,b,l/b) , 



;i2) 



where 

P(x, 6, 1/6) = <P{x, h) + 0K(l/6)) . (13) 

The evolution of V in 6, denoted by 0(as), has been neglected. The h depen- 
dence in 0, corresponding to the intrinsic transverse momentum dependence 
of the pion wave function [^, ^ , will not be considered here. 






Similarly, the renormalization group analysis applied to H gives 

H{xi,b,Q,j2) = exp 

xH{x„b,Q,t), (14) 

where t is the largest mass scale involved in the hard scattering, 

t = m.ax{y/xiX2Q, 1/6) . (15) 

The scale ^xxX^Q is associated with the longitudinal momentum of the hard 
gluon and 1/6 with its transverse momentum. Combining all the exponents 
derived above, we obtain the lowest-order expression for the pion form factor, 

/■l /"OO 

^AQ^) = ISttCf / dxidx2(j){xi)(f){x2) bdbas{t) KQ{JxiX2Qh) 
Jo Jo 

X exp[-S'(a;i, X2, 6, Q)] , (16) 

where the complete Sudakov logarithms are given by 

2 2 t 

S(xi, X2, b,Q) =Y, [s{xi, 6, Q) + s(l - Xi, 6, Q)] - — In^ , (17) 

i=i Pi 6 

with i = ln(t/A). Cp is the color factor defined by tr (T"T") = NcCp, Nc 
being the number of colors. Kq is the modified Bessel function of order zero, 
which is the Fourier transform of the gluon propagator to the 6 space. 

Variation of S with 6 has been displayed in |T[], which shows a strong 
falloff in the large 6 region, and vanishes for 6 > 1/A. Hence, Sudakov 



suppression selects components of the pion wave functions with small spatial 
extent b, and makes the hard scattering more perturbative. If b is small, 
as with its argument set to t as in eq. ([T5|) will be small, regardless of the 
values of the a;'s. When b is large and xiX2Q'^ is small, as is still large. 
However, the Sudakov factor in eq. (|16D strongly suppresses this region. Since 
the main contributions to the factorization formula come from the small b, 
or short- distance, region, the perturbation theory becomes relatively self- 
consistent. It is also easy to show that the modified factorization formula 
including Sudakov corrections reduces to the standard one as 6 ^ 0, where 
the important logarithms diminish. 



3. The Decay B -^ nliy 

We now generalize the above analysis of the pion form factor to the ex- 
clusive processes involving both light and heavy mesons. In particular, we 
work with the semi-leptonic decay B -^ ttIu. PQCD is appropriate to this 
process when the pion is energetic, because the b quark mass provides a 
large scale. We shall analyze the leading regions of radiative corrections to 
such a process, and derive the Sudakov factor including both leading and 
next-to-leading logarithms. The amplitude of this decay is written as 

A{P,,P2) = ^V^,u^,{l - ^,)l{7r{P2)\bru\B{P,)) , (18) 

where the four-fermion interaction with the Fermi coupling constant Gp has 
been inserted. Pi and P2 are momenta of the B meson and the pion, respec- 
tively. We start with the lowest-order factorization for the matrix element 
M^ = {7T{P2)\b'~f'^u\B{Pi)) with an exchanged hard gluon as shown in fig. 3, 
the left-hand side being the B meson at rest and the right-hand side a fast- 
recoiling pion. The heavy b quark is represented by a bold line. The symbol 

10 



X denotes the electroweak vertex with the KM matrix element Vub, from 
which a lepton pair emerges. 

Parton momenta are assigned as in fig. 3. The b quark carries Pi — ki, and 
the accompanying hght quark carries ki. They satisfy the on-shell conditions 
(Pi — kiY ^ ml, Pf = m\ and /c^ ~ 0, ms being the B meson mass. In 
the Breit frame Pi has the components P^ = Pf = rriB/V^ and vanishing 
transverse components, ki may have a large minus component with small 
amount of transverse components ki^, which will serve as the infrared cutoff 
of Sudakov corrections below. The assignment of parton momenta on the 
pion side is similar to that in the case of the pion form factor as shown in 
fig. 3. The large component of P2 is P.^ = r]mB/V^, rj being related to 
the energy fraction of the pion by P2 = rim/2. The physical range of r] is 
< r] < 1, since the pion carries away at most half of the rest energy of the 
B meson. The transverse momentum associated with the valence quarks of 
the pion is denoted by Vi2T- The invariant mass of the lepton- neutrino pair 
produced in this decay is given by ml = (Pi — P2)^ = (1 — ffjm^. 

We then consider radiative corrections to the basic factorization for the 
heavy-to-light transition. The essential step is again to find out the leading 
regions of radiative corrections in axial gauge. For reducible corrections on 
the pion side, the conclusion is the same as before: they produce double 
logarithms with soft ones cancelled in the asymptotic region, and can be 
absorbed into the pion wave function, which give rise to the evolution of 
the wave function. Irreducible corrections, with an extra gluon connecting a 
quark in the pion and a quark in the B meson, give only soft divergences, 
which also cancel asymptotically. Hence, they are absorbed into the hard 
scattering amphtude. 

On the left-hand side, three diagrams showing the 0{as) reducible ra- 
diative corrections, are exhibited in fig. 4. Fig. 4a, giving the self-energy 
correction to the massive h quark, produces only soft enhancements, and is 
thus not leading. If /cf is small, collinear divergences in figs. 4b and 4c, which 
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arise from the loop momentum with a large component parallel to ki, will 
not be pinched, and they also give only soft enhancements. This is consistent 
with the physical picture that the soft light quark can not interact with the 
heavy quark through a fast moving gluon. However, from the B meson wave 
functions given in section 4, it is observed that there is substantial proba- 
bility of finding the light quark with fcf of order m^, even though the wave 
functions peak at small ki . Therefore, figs. 4b and 4c contribute collinear 
enhancements. Note that fig. 4b gives soft divergences which are not com- 
pletely cancelled by those from figs. 4a and 4c even in the asymptotic region. 
In conclusion, figs. 4b and 4c indeed contain double logarithms, which must 
be organized by the resummation technique. 

Since the collinear enhancements on the B meson side are less impor- 
tant due to the suppression of the wave functions, reducible corrections are 
basically dominated by soft enhancements, and can be absorbed into the B 
meson wave function, which is also dominated by soft dynamics. This ab- 
sorption of reducible corrections is similar to that on the pion side, where 
reducible corrections are dominated by collinear enhancements. Therefore, 
the factorization picture in fig. 3 still holds after radiative corrections are 
included. 

With the above reasoning, we can write down the factorization formula 
for the decay B -^ nlu in the transverse configuration space, 

M^ = X dx,dx,y^^P.(x2,b„P„/.) 

X H^{xi,X2, bi, b2, m, fi) Vb{xi, bi. Pi, /i) , (19) 

where both the pion and B meson wave functions, Vt^ and Vb, contain leading 
double logarithms. Here we have to introduce two 6's, because the virtual 
quark line in the hard scattering may not be far from mass shell, and can not 
be shrunk to a point. This is contrary to the case of the pion form factor, 
and detailed explanation will be given later. Therefore, we need bi to denote 
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the separation between the two valence quarks of the B meson, and 62 for the 
pion. The approximation m^ ^ rriB = ttl has been made. The momentum 
fraction Xi is defined as k^ / P{ . H^ is the Fourier transform of the hard 
scattering amphtude derived from fig. 3, whose exphcit expression will be 
given in section 4. The resummation of the double logarithms in P^ has 
been performed in the previous section. Here we quote the results directly 



with Q set to rjm and ri/ set to 4 |]T3|. In the below we shall concentrate on 

Vb. 

There are two major difiiculities in summing up the double logarithms in 
figs. 4b and 4c. Firstly, fig. 4a, which gives only single logarithms, must be 
excluded. Secondly, there are many invariants that can be constructed from 
Pi, ki and n such as P^, Pi ■ ki, Pi ■ n, ki ■ n and n^, which are involved in 
Vb- In the pion case the number of invariants is much less. There are only 
P-n and n^, and the resummation is thus simpler. The fact that Vb contains 
many invariants fails the technique of replacing d/dm by d/dn as introduced 
in section 3, because in this case some large scales like Pf can not be related 
to n. 

However, the above difficulities can be overcomed by applying the eikonal 
approximation to the heavy quark line as shown in fig. 5. In the coUinear 
region with the loop momentum parallel to ki and in the soft region, the b 
quark line can be replaced by an eikonal line: 

(/^i-A+^ + m)7" ^^ (Pi-fci)" . . 

(Pi -ki + q^-m^ ^ {Pi-ki)-q ' ^ ^ 

where the remaining part R either vanishes as contracted with the matrix 
structure of the B meson wave function, or is less leading. The factor l/[(Pi — 
ki) ■ q] is associated with the eikonal propagator, and the numerator (Pi — 
ki)" is absorbed into the vertex, where a gluon attaches to the eikonal line. 
The physics involved in this approximation is that a soft gluon or a gluon 
moving parallel to ki can not explore the details of the b quark, and its 
dynamics can be factorized. The idea is similar to that employed in HQET. 

13 



An explicit evaluation of radiative corrections confirms this approximation. 
By this means, the first difficulty is resolved, because self-energy diagrams 



of an eikonal line are excluded by definition |T6 . 

The eikonal approximation also reduces the number of large invariants 
involved in Vb- We have the scale invariance in Pi — ki ^ Pi as shown by 
Feynman rules for an eikonal line in eq. (|20| ) , which corresponds to the flavor 
symmetry in HQET, in addition to the scale invariance in n. Hence, Pi does 
not lead to a large scale, and the only remaining large scale is fcf , which must 
appear in the ratios (fci ■ nY /ri^ and {ki ■ PiY / P^. An explicit lowest-order 
investigation shows that the second scale {ki ■ PiY / Pf in fact does not exist. 
Therefore, with the eikonal approximation the problem is simplified to the 
one in analogy with the light-meson case. Now Vb depends only on the single 
large scale z/'^ = {ki ■ nf jr? ^ and d/d/cf can be replaced by d/dn. 

Following the same procedures as in section 2, the differential equation 
of Vb is derived as 



-Vn = -T^^V 



B 



IC{bf,) + ^Qiu'/f,) 



Vb , (21) 



d In 771 d In fci 

where the lowest-order /C can be obtained from fig. 6a, and Q from fig. 6b, 
with the square vertex representing 

77 
gT- fc^ . 

ki ■ nq ■ n 

Note the absence of the diagrams corresponding to self-energy corrections to 
the eikonal line. 

Comparing fig. 6a with 2a, we find that the evaluation of /C for the B 
meson is similar to that for the pion except the third diagram. This extra 
diagram is finite without ultraviolet and infrared divergences, as justified by 
its integral 

^^^J (Anrki-nq-nPi-q q^ ^ ' ^ ^ 
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Hence, it does not spoil the renormalizat ion-group invariance of the Sudakov 
logarithms. For the specific choice ra oc Pi as in the case of the pion form 



factor, it is easy to show that (p^) vanishes. Therefore, the functions fC and 
Q for the B meson are in fact the same as those for the pion. 

Substituting the expressions of /C and Q into eq. (^1]), we obtain the 
solution 



Vb{xi, 6i, Pi, /x) = exp [-s{xi, 61, m)] Vb{xi, h, /x) 



(23) 



where the exponent s is given by eq. (|^) but with nf = 4 |]T^. Summing up 
single logarithms in VB{xi,bi, fi), eq. (p3D becomes 



'Pij(xi,6i,Pi,/i) = exp 



-s(xi, hi,m) 



2/' ^-^iMm 
Ji/bi fi 



X(f)B{xi,bi) + 0{asil/bi)) 



(24) 



where the anomalous dimension 7^ is the same as before. Again, the evolution 
of Vb in 61 and the intrinsic 61 dependence of the wave function (pB will be 
neglected below. Including the summation of single logarithms in iP^ and 
the results from P^, we obtain the complete Sudakov exponent 



S{xi,bi,r],m) = s{xi,bi,m) + s{x2,b2,r]m) + s{l - X2,b2,r]m) 

-— (in^ + ln^ I . 
A V h b2J 



(25) 



The variable 6j is defined as in eq. (|^), and t is the largest scale involved in 
the hard scattering, which will be discussed in section 4. Similar expression 
to eq. ( p5D for the pion form factor, which also involves two 6's, has been 
obtained in HI. 



4. Numerical Results 
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Having derived the Sudakov factor for the process B -^ tt/z/, we now 
evaluate its decay rate, and examine how much contribution comes from the 
perturbative region with small 6j. The explicit formula for lowest-order if^^*^) 
from fig. 3a is written as 



^li{a) 



tr 



X 



la 



75 /P2_u /^i ~ ^2 /P2+ A2T + m ^ (/Pi + m)75 



=7- 



V^K' {Pi 



X2P2 + ^2tY - m? 



V2K 



{X2P2 -h+ k2T)2 

4(1 + X2'r])g'^CF'm^ 
[x2rim'^ + \<-2t] [xiX2f]'m'^ + (ki^ — k; 



2TJ 






(26) 



where the factors 75 /p2/V2A^c and (/Pi + 771)75/ v^2iVc come from the matrix 
structures of the pion and B meson wave functions, respectively |TB[. The 
relation k^ = xim/\/2 has been inserted. In the second expression the k2T 
dependence in the fermion propagator is not neglected. For the pion form 
factor the corresponding transverse momentum dependence is negligible, be- 
cause there is a factor X2 in the numerator, which cancels the singularity 
from X2 — i> in the denominator. However, for the B meson decays, due to 
the massiveness of the b quark, such a cancellation does not appear as shown 
in eq. (PB]). To ensure that the virtual quark be part of the hard scattering, 
k2T must be retained. Therefore, with the inclusion of transverse momenta, 
we need not to perform the substraction of an on-shell fermion propagator 



from the hard scattering as in [13 



Similarly, the expression for lowest-order H^^^^ is given by 



fftJ-ib) 



tr 



75 /P2 

7a /7nr^ 7 



/'2+ A ^ J/Pi + m)75 



r7 



X 



Vm' {P2 + h)- 

{X2P2 -h + k2T)2 



V2K 



[xirjm? + kfj,] [xiX2rim? + (ki^ - k2T) 
Ag'^CpXim^ 



-Pf 



[xirirn? + k{j][xiX2rim? + (kiT - k2T) 



rPo^ 



(27) 
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To derive the second expression, we have replaced ki by 

Here the transverse momentum ki^ in the fermion propagator is neghgible, 
because the singularity from xi ^ is removed by the numerator. However, 
we will keep it for consistency. 

Performing the Fourier transform of eqs. (|26| ) and (|27|) , the matrix element 
M'^ can be written as 

M^ = hPr + f2P^ , (29) 

in which the factorization formulas for the form factors /i and /2 are given 
by 

/•l roc 

/i = lOvrCj^m^ / dxidx2 / &id6i62d62 0B(2;i)07r(a;2) 

xxir]h{xi,X2, bi, 62, Vj "^) exp[-S'(xi, 6^, rj, m)] , (30) 

and 

rl /"OC 

/2 = IGrcCFm^ / da;ida;2 / bidbib2db2 (j)B{xi)(j)Tr{x2) 
Jo Jo 

X [-Xih{xi,X2, 61, 62, ??, m) + (1 + X2r])h{x2, Xi, 62, ^i, ??, "i-)] 

X exp[-S'(a;j, 6i, r], m)] , (31) 

respectively, with 

h{xi,X2,bi,b2,i],'m) = as{t)Ko{y/xiX2i]mb2) 

x[e{bi - b2)KQ{.JxYnmbi)h{^/x^mb2) 

+e{b2 - 6i)iro(v/iI77m62)/o(v/^^^i)] -(32) 



Jo is the modified Bessel function of order zero. From eqs. (^61) and (|27|) , we 
choose the largest scale t associated with the hard gluon as 

t = max (y^x 1X2 77m, 1/61, 1/&2) • (33) 
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Basically, the above expressions for /« are similar to that for the pion form 
factor but with different Sudakov logarithms. 

To evaluate fi, we consider the following two models of (psi^), which have 
been adopted in [0. They are the oscillator wave function ||l] 



NJx{1 — x) exp 



m. 



2cj2 I 2 2m2 



'■B , 



NJx{1 — x) exp 



ni 



2cj2 



-X 



in our approximation rri^ ?« m^, and [|I 



$ 



W, 



B 



x,kT) 



C + 



mt 



+ 



JV'T"' 



n -2 



(34) 



(35) 



1 — X ' x(l — x) 

where x is the momentum fraction of the light quark in the B meson. The 
parameters in 0g are u = 0.4 GeV and m = 5.28 GeV. The constant N is 
determined by the normalization of the wave function: 



dn, 



h 



dx(/)V(a;) - „ /Tj , 
2v3 



(36) 



/b = 160 MeV being the B meson decay constant [T^, which leads to A^ = 
2.24 GeV. 

Similarly, the constants A^' and C in $^ are determined by the following 



normalizations 18 



/ dx f d^kT<l>^i^\x,kT) 



I 



B 



2v^' 
'dxJd'kr[<^^i'\x,kT)r = ^, (37) 

1.232 GeV^ and C = -0.99993771^ are obtained. The 

{II) . 
B gives 

<P^i'\x,b) = J d'kT<^^i'\x,kT)e^'^--'' 
nN'bx^(l~x? 



from which A^' 
Fourier transform of $ 



Jm?x + Cx{ 



'^ — Ki Uw?x + Cx{l - x)h] , (38) 
1 — x) ^ ^ 
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with Ki the modified Bessel function of order one. As stated before, we 
neglect the intrinsic b dependence of the wave function, and define 

4'"(.)^ii;j.r(M)^£|^:^^. (39) 

Obviously, both of the models peak at small x, which signifies the soft dy- 
namics in the B meson. However, the probability at intermediate x is indeed 
comparable at least in model II, and the resummation of double logarithms 
performed in section 3 is essential. At last, we adopt the Chernyak-Zhitnitsky 
model |2^ for the pion wave function: 

Mx) = 5V3Ux{l-x){l-2xy (40) 

with /^ = 93 MeV the pion decay constant. 

We are now ready to evaluate /i and /2 numerically in eqs. ( ppp and (PT|). 
Results of fi for the two models of (pB, with bi and 62 integrated up to the 
same cutoff bic = 62c = be, are shown in fig. 7. It is found that at r/ = 0.3 
approximately 50% of the contribution to fi comes from the region where 
as{l/bc) < 1, or equivalently, be < 0.5A. At 77 = 0.4, 55% of the contribution 
is accumulated in this perturbative region. As 77 = 1, perturbative contri- 
bution has reached 70%. It implies that the PQCD analysis of the decay 
B -^ ttIu in the range of 77 > 0.3 is relatively reliable according to the crite- 
ria given in [Q, Q . Therefore, we shall accept the modified PQCD predictions 
for 1] > 0.3, and use them to evaluate the differential decay rate. It is also 
observed that the expressions with 0^ employed are more perturbative, 
because the large /cf region is less suppressed by the wave function, and the 
double logarithms are stronger. The predictions from 0^ are larger, because 
0^ enhances the contribution from the end-point region of Xi. 



In the approach of |]T3[ where the transverse momentum dependence was 
not included, instead, energies of the virtual quark and gluon involved in the 
hard scattering were taken as the ultraviolate and infrared cutoffs of radia- 
tive corrections, respectively. The resulting Sudakov logarithm proportional 
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to Inm/ki, however, gives weaker suppression. From the steepest descent 
approximation of their Sudakov factor, the saddle point at fcf around 1.4A 
for r/ = 1 has been found. Using the same method, we determine the saddle 
point of our Sudakov factor in eq. (|25D by 

dS dS 

from which a larger scale 1/bi = lOA for A = 0.1 GeV, or 1/bi = 6 A for A = 
0.2 GeV is obtained. It is then obvious that the perturbative expansion in our 
approach is more reliable. At this large scale the double logarithms contained 
in the 0{as) radiative corrections to the hard scattering with a triple gluon 
vertex, which have been considered in |T^, are in fact not important. 



We observe that the magnitude of /2 is much larger than that of /i, 
especially in the small rj region. This fact is consistent with their behaviors 
in the soft pion limit as obtained in the framework of HQET ||12|, in which 
/2 is found to have a pole at r/ — ;► 0: 

lim/2^^^BB.. , (42) 

for rriB* ~ ruB, ruB* being the B* meson mass. In the above expression 
Qbb'tt is the BB*n coupling constant, and /b* the decay constant of the 
B* meson. Assuming t] = 0.3, /b* ~ 1-1/b and gBB'Tv ~ 0.75 |22 



m 



eq. (|4^), we obtain /2 = 9.5 which is close to our predictions at r^ = 0.3. 
Certainly, this extrapolation of the soft pion theorems may not be reliable, 
but it is interesting to observe that it is consistent with the PQCD results at 



the middle value of rj. /i in the soft pion limit vanishes like 1 — JniB* /niB 

m- 

With results of /i and /2, we can compute the differential decay rate of 
B^ -^ Ti^l^u with vanishing lepton masses [|T^: 

dr _ 2G|mV|. , .|2 

ar] YooTT'^ 

= \Vub\^R{v), (43) 
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for 7] > 0.3, where the second formula defines the function R{t]). Results 
of R{r]) are shown in fig. 8, from which the decrease of dT/dr] with rj is 
observed. The behavior predicted here is opposite to that given in [|1^ 



which shows an increase in t] starting with zero at rj = 0.5. Such a dip at 
the middle value of i] is attributed to the substraction of an on-shell fermion 
propagator from the hard scattering |jl3|, which is, however, not necessary in 



our treatment. Furthermore, the predictions in |T^ for the differential decay 
rate are almost 10^ times smaller than ours. The reason is again traced back 
to the substraction of the fermion propagator. Therefore, the transverse 
momentum dependence plays an essential role in the analysis of B meson 
decays. 

The differential decay rate in the soft pion limit can be obtained using 
eq. (H): 

lim R{rj) = ^ha^lkgl , (44) 

which shows a linear relation with tj. Extrapolating eq. (|^) to intermediate 
7] as shown in fig. 8, we observe a fair match between the soft pion and 
PQCD predictions around rj = 0.3. It implies that the modified PQCD 
formalism is successful at large rj, but becomes worse quickly in the soft 
pion limit. Similarly, the soft pion technique is appropriate at small rj, but 
gives an overestimation to the considered process in the perturbative region. 
The overlap of these two approaches indicates the transition of the B meson 
decays to PQCD at middle values of rj, and the complementarity between 
the soft pion theorems and modified perturbative formalism [^ ^ . 

We then estimate the total decay rate F by integrating dT/dr] using 



eq. (^) for rj < 0.3 and using our PQCD predictions for rj > 0.3. We obtain 
0.6 X lO^^^lKfeP GeV from the soft pion theorems, and 1.9 x 10"^^|KfeP and 
0.7 X lO^^-'^lVubp GeV for the use of 0^ and 0^ , respectively, from the 
modified PQCD formalism. Their sum gives F ^ 2.5 x 10~^^|\4fep GeV for 
model I, and 1.3 x 10~^^|\4bp for model II. They correspond to branching 
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ratios 0.5 x lO^lKbP and 0.26 x lO^lKbP, respectively, for the total width of 
the B^ meson is (0.51 ± 0.02) x 10^^ MeV [Q. Current experimental limit 



on the branching ratio of 5° ^ vr l^u is 3.3 x 10 ^ ||2^. We then extract the 



the matrix element \Vub\ < 2.6 x 10 ^ from model I and \Vub\ < 3.6 x 10 ^ 
from model II. These upper limits are close to the value 0.003 given in the 



literature 24 



5. Conclusion 

In this paper we have applied the resummation technique to the semi- 
leptonic decay B -^ vr/z/, and derived the Sudakov factor up to next-to- 
leading logarithms in this heavy-to-light process. The idea is to apply the 
eikonal approximation to the heavy b quark line such that its nonleading self- 
energy diagram is excluded, and the number of large scales involved in the B 
wave function is reduced. The resummation of double logarithms in the heavy 
meson is then simplified to the one in analogy with the light meson. The 
modified PQCD calculation of the differential decay rate including Sudakov 
effects has been examined and found to be reliable for t] above 0.3. By 
combining our predictions with the soft pion results and comparing them 
with experimental data, we have estimated the total decay rate, from which 
the upper limit 0.003 for the KM matrix element \Vub\ is obtained. 

We do not observe the dip at ?7 = 0.5 for the differential decay rate as 



predicted in |]13|, which arises from the substraction of an on-shell fermion 
propagator from the hard scattering. This substraction is not necessary in 
our analysis because of the inclusion the transverse momentum dependence. 
The behavior in rj is similar to that of the pion form factor in Q^, which 
is opposite to the results of [|r^. We do not consider the evolution in b 
and the intrinsic transverse momentum dependence of wave functions in this 
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work. However, we estimate that these two modifications cancel at least 
partially, since the former gives an enhancement [Q to the decay rate, but the 
latter leads to a suppression |jl5[. Certainly, this subject needs more detailed 



investigation. Our formalism can be easily applied to a similar semi-leptonic 



decay B -^ plu |jT^ and other non-leptonic B meson decays, which will be 



published elsewhere. 
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Figure Captions 

Fig. 1. 0{as) radiative corrections to the basic factorization of the pion 
form factor. 

Fig. 2. Graphic representation of eq. (1). 

Fig. 3. Lowest-order factorization for the decay B -^ tiIu. 

Fig. 4. 0{as) radiative corrections to the B meson wave function. 

Fig. 5. Eikonal approximation for the b quark hne. 

Fig. 6. Lowest-order diagrams for (a) the function JC and (b) the function 
Q associated with the B meson. 

Fig. 7. Dependence of (a) /i and (b) /2 on the cutoff be derived from (p^ 
(sohd hues) and from 0^ (dashed hues) for (1) r] = 0.3, (2) r] = 0.4, and 
(3) r] = 1.0. 

Fig. 8. Dependence of R{r]) on i] derived from (1) 0^ and from (2) (p^ 
(sohd hues). Results from the soft pion theorems are also shown (dashed 
line) . 
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